A 2 PN/RM metric of General Relativity 
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We present a derivation of 2PN/RM metric field equations from the Einstein field equation in 
General Relativity. We use the exponential parametrization and the c~ A isotropic spatial coordinates 
such as in IAU2000 recommendations. 
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I. INTRODUCTION 

In a foreseeable future, telemetry will reach a new ac- 
curacy using laser links and high accuracy clocks (such 
as the TIPO experiments which is in study 1 ). At 
such experiments accuracies, the effective 1,5 PN/RM 
(acronym defined in section TTb metric defined from the 
1PN IAU2000 [2] metric will be no longer enough to fit 
data (2 PN/RM effects lead to terms of about KT 10 s for 
a flying time of about 1000s when tagging units accuracy 
reach at least 10~ n s). Hence, it is necessary to develop 
the metric up to the 2 PN/RM order. Our work is com- 
pletely based on General Relativity (GR) and hence we 
don't consider any alternative theory in this paper. This 
means that our result would have a practical usefulness 
in telemetry only if the GR theory is the "true" theory of 
gravity or if alternative theory deviates from GR only by 
sufficiently small post-Newtonian numerical parameters 
in the metric. This is compatible with present available 
data. 

We start our work from the "exponential parametriza- 
tion" such as in DSX |3] or in the IAU2000 resolutions 
[2]. Hence, our notations come from [3]. Other works 
on the subject have already be done ([4], [5]) but they 
didn't write their metric in the convenient exponential 
parametrization and these works have been done in har- 
monic gauge from the start to the end. This is not our 
case since we fix the gauge only at the end when we want 
peculiar solutions in some peculiar gauges. 
In section [H] we introduce a new terminology before re- 
calling some well known results from the DSX paper 
[3] in section III Then we derive Einstein's equations 



in our metric's parametrization and finally we give for- 
mal solutions in any gauge that respect the exponential 
parametrization - with a particular attention to the har- 
monic gauge. 



II. TERMINOLOGY : DEFINITION OF THE 

PN/BM AND PN/RM METRICS 



order \jGMjr or less, M being some caracteristic mass 
of the system) for both the sources and the (test) body. It 
formally consists in looking for solutions under the form 
of an expansion in powers of 1/c. The usually so-called 
nPN order terms in the metric, leading to c~ 2n terms in 
the equation of motion of a body describing a bounded 
orbit, are terms of orders c _2n ~ 2 in g Q0 , c~ 2n ~ x in g oi 
and c _2n in g^. In this paper, a metric developped this 
way will be refered as the nPN/BM metric (BM meaning 
"Bounded Motion" for test particles). It is particularly 
well-adapted for studying bounded motions in systems 
made by non-relativistic massive bodies, as the Solar Sys- 
tem is. 

However, since we are interested in the propagation of 
light, we are lead to relax the hypothesis on the veloc- 
ity of the test particle whose motion is considered. Of 
course, this doesn't change the full metric, but the terms 
to be considered in the metric components are not the 
same as in the PN/BM problem. Indeed, the terms lead- 
ing to c -2 " terms in the equation of motion of a test 
particle moving with relativistic velocity are terms of or- 
der c~ 2n in g a p (ie. in both g m , g oi and g^). In this 
paper, a metric developped this way will be refered as 
the nPN/RM metric (RM meaning "Relativistic Motion" 
for test particles) . A PN/RM metric is particularly well- 
suited for studying relativistic motions of test bodies (for 
instance, the propagation of light) in systems made by 
non-relativistic bodies, as the Solar System is. 

PN/BM and PN/RM orders are illustrated in figure [T] 
The present paper deals with the PN/RM problem 
since we are concerned in propagation of light. 



III. 1,5 PN/BM METRIC 

We recall here results from DSX^ on the 1,5PN/BM 
metric. First, the "exponential parametrization" with 
the spatial isotropy condition makes the metric writes 



The PN approximation is based on the assumption of 
a weak gravitational field and weak velocities (ie. of the 
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IV. 2 PN/RM METRIC 
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At the 2PN/RM order, we get 6 new degrees of free- 
dom in the metric. As we know, from the 10 degrees of 
freedom in the equation of GR, there are only 6 cleared 
by the 10 equations since 4 degrees of freedom correspond 
to a gauge invariance. Indeed, the metric is determined 
by the field equations modulo the knowing of the used 
coordinates. Since, there are four coordinates, there are 
4 degrees of freedom not reduced by the 10 field equa- 
tions. As seen in the previous section, to the 1,5PN/BM 
approximation correspond 4 metric field equations and 
the gauge invariance, characterized by A, keep 1 degrees 
of freedom. Then, we expect from the 2PN/RM order 6 
equations with a gauge invariant field (characterized by 
a 3- vector) which (the equations) reduce the 6 degrees of 
freedom to 3. Let us fix the covariant form of the metric 



FIG. 1: General scheme of orders taken into account in 
PN/BM and PN/RM metrics. 



With jij = Sij + 0(c 4 ). From which follows 
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The Einstein field equation then reduces to four equa- 
tions on w and w. 
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.dtdjWj = -4irGa + O^) (7) 
i - d t diw = -4ttG<7' + 0(c- 2 ). (8) 
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a 1 = c- x T Qi and A = 5-8. 



Hence, to the 1,5PN/BM approximation correspond 4 
equations of the GR equation and the gauge invariance 
keep 1 degree of freedom (indeed <r'' is approximately 
divergence free). This gauge invariance is characterized 
by the arbitrary differentiable function A in DSX. 
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The IAU2000 resolution then corresponds to a peculiar 
choice of gauge in order to fix the last degree of free- 
dom : the harmonic gauge. This condition writes, in the 
1,5PN/BM approximation {g a0 T° ap = 0(c- 5 ),g a ^T l a/3 = 
0(c~ 4 )}. The four metric field equations then reduce to 





D m w = 


-iirGa + 0(c- 4 ) 


(11) 




Aw, = 


-A-kGct 1 + 0(c- 2 ) 


(12) 


where D m 


= a - a 2 . 
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Then we can derive the 2PN/RM field equations from the 
Einstein field equation R a P = ^f- (T a P - \g al3 g P aT pa ). 



-j@ij(T k i) + Zij(w,w k ) = 

8nG (IrfJT 00 4- T^ - -n^T kk - — 
c 4 \2 V 2 V c 2 



jjjoo ( 19 ) 



Where 



% M = d 2 ki T kj + d 2 kjTik - A nj - d 2 T kk , (20) 



Eij (w, w k ) = 

-^n m w+— (d^Wj + dojWi) — -^diwdjW-] — ^-wAw. 

(21) 

Then, considering w and w k as source terms and iden- 
tifying the relation -^ (^-wAw) = ^f-r/^wT 00 + 
0(c~ 2 ) from (11), we get 



9« (m) = 47rG(T ij ' 



[(9 t j w j + ^w, — diwdjw) 



2rf* (d 2 
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Where ct" = T % ' 3 — °rfJT kk . As expected, the oper- 
ator is gauge invariant through a relation on any 3- 
vector field. Indeed, we have 0jj (tui + d k Ai + diAk) = 
Q»j ( r fcO: whatever is the 3- vector field A&. Hence, the 
2PN/RM gauge invariance is characterized by both the 
arbitrary differentiable function A and the arbitrary dif- 
ferentiable 3- vector A k . 
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T 'ki = T ki + d k Ai + d[A k + -S k id t X. 
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level, the harmonic constraint is g a ^^' 1 a a = 0(c 5 ). That 
gives the harmonic conditions on the "potentials" at the 
2PN/RM level 



d t w + d k w k = 0(c 2 ) 

d k T lk -\ dl Tkk + d tWl = 0{c^). 
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Then 0y (Tjy) writes: 



0y (m) 



-{A Tij + OtidiWj + djWi)}. (30) 



Hence, considering w and Wi as source terms, the solu- 
tion of the Einstein field equations in the harmonic gauge 
is : 



A" 1 {-4ttG(t ij - diwdjw} 



(31) 



A. General formal solution 

We can show that it always exists a gauge transforma- 
tion r kl — Tki + d k Ai + diA k that reduces the operator 
6 to a Laplacian (0^- (rjy) = —AtL). To be placed on 
this gauge, A k must satisfy AAj — d k r k j — \djT kk which 
is clearly inversible since we know the laplacian's green 
function. Hence, there always exists a peculiar gauge 
where the solution is 

th = A _1 {-47rG<7 4 * 

+ [(dlw, + %jWi - d iW d jW ) + 2tj« (d 2 w + d 2 tk w k )] }, 



where 



A- 1 {f(t,x)} = 
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Then, from this equation, we can get the solution T k \ in 
any gauge, using any vector field V k with r kl = tu + 
dkVt+dM. 



B. The harmonic gauge 



The harmonic gauge condition stands g^T^g = 
(Or d a {\f—gg ai ) = with g the metric determinant). 
As seen in DSX [3], the harmonic 1,5PN/BM condition 
9 a ^^aB = 0(c~ 4 ) gives the usual harmonic constraint on 
the "potentials" (d t w + d k w k = 0(c~ 2 )). While, the con- 
ditions g a ^ % a a — 0(c) give no more information due 
to the "isotropy condition" which makes spatial coordi- 
nates always harmonic modulo 0(c~ 4 ). At the 2PN/RM 



where us and uij satisfy the 1,5PN/BM harmonic condi- 
tion d28b. 



from non-harmonic to harmonic gauge 

If Tij is not an harmonic solution, then the gauge 
transformation t,-, 



harmonic solution if Aj satisfies 



e ij + d{Aj + djAi makes ey an 



AAj = d k T kj - -djTkk + d t w.j 



(32) 



where Wj must be expressed in the harmonic gauge with 
the 1,5PN/BM field equations (12). Since it is inversible, 
this relation always allow us to put the metric in the 
harmonic gauge, whatever is the gauge from which we 

start. 



V. CONCLUSIONS 

In order to fit future inter-planetary telemetry data, we 
developed a GR metric up to the 2 PN/RM order with 
the DSX "exponential parametrization" . We gave both 
general and harmonic formal solutions. A more detailed 
version of this result can be found in [S] . An explicit solu- 
tion of these equations and of the (non-massive-)Scalar- 
Tensor 2PN/RM field equations in the case of the solar 
system is in preparation. 
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